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ABSTRACT 


This thesis provides a comparison of the accuracies of two methods for fitting a 
Weibull distribution to a set of aircraft engines time-between-failure data. 

One method used is the Maximum Likelihood Method and assumes that these engine 
failure times are independent. The other method 1s a Modified Method of Moments procedure 
and uses the fact that if time to failure T has a Weibull distribution with scale parameter 4 and 
shape parameter B, then T? has an exponential distribution with scale parameter 4°. The latter 
method makes no assumption about independent failure times. 

A comparison is made from times that are randomly generated with a program. The 
program generates times in a manner that resembles the way in which engine failures occur in 
the real world for an engine with three subsystems. These generated operating times between 
failures for the same engine are not statistically independent. This comparison was extended to 
real data. 

Although the two methods gave good fits, the Maximum Likelihood Method produced 
a better fit than the Modified Method of Moments. Explanations for this fact are analyzed and 


presented in the conclusions. 


THESIS DISCLAIMER 
The reader 1s cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While effort has been made, within the time 
available, to ensure that the programs are free of computational and logic errors, they cannot 
be considered validated. Any application of these programs without additional verification is at 


the risk of the user. 
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EXECUTIVE SUMMARY 


This thesis provides a comparison of the accuracies of two methods for fitting a 
Weibull Distribution to a set of aircraft engines time-between-failure data. 

One method used is the Maximum Likelihood Method and assumes that these engine 
failure times are independent. The other method is a Modified Method of Moments procedure 
and uses the fact that if time to failure T has a Weibull Distribution with scale parameter 4 and 
shape parameter B, then T° has an exponential distribution with scale parameter 2°. The latter 
method makes no assumption about independent failure times. 

In a fleet of fielded aircraft engines, engines are repaired when a fault is detected. In 
addition, some critical components are inspected/repaired after a designated number of flight 
hours. When faults are detected the engine is repaired and the component is replaced on the 
engine unless the engine has accumulated a total number of flight hours near to the hard 
inspection time of a critical component(s). If this is the case the critical components(s) are also 
inspected/repaired. In nearly all cases, the engines are not completely overhauled to restore 
them to new condition. After repair, components in the engine have mixed ages with respect to 
flight hours with some components new or repaired. Consequently the times between failures 
on a specific engine are not independent. 

A comparison is made for times that are randomly generated from a program. The 
program generates times in a manner that resembles the way in which engine failures occur in 
the real world for an engine with three subsystems. These operating times between failures for 
the same engine are generated in a way that they are not statistically independent. This is 
accomplished by using a program in S-Plus which simulates failure times for an engine with 
three subsystems. The failure time distributions of each one of these subsystems are Weibull 
with fixed parameters determined a-priori and based on real world experience. The simulation 
gathers three failure times for each engine. The study includes failures for two different sets of 


20 and 100 engines respectively. Once this comparison was completed, the study was extended 


to real data using TF-34 engine operating times between failures from the Naval Aviation 
Logistics Data Analysis (NALDA) data base. 

The estimation of the parameters for the Weibull Distribution which fits the failure 
times was done using spread-sheet software programs. The accuracy of each method was 
measured using the Kolmogorov-Smimov test statistic, embedded in the S-Plus software 
package. 

Although the two methods gave good fits, the Maximum Likelihood Method produced 
a better fit more often than the Modified Method of Moments. Explanations for this fact are 


analyzed and presented in the conclusions. 
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I. INTRODUCTION 


A. BACKGROUND 


ie Engine Reliability Analysis Program (ERAP) 

The Naval Air Warfare Center Aircraft Division Propulsion and Power Engineering 
(NAWCADPPE) was tasked with analyzing aircraft engine failure data to estimate aircraft 
engine reliability and maintainability factors and to assess the effectiveness of existing naval 
aircraft engine maintenance policies. The ERAP program was established to help accomplish 
this task. It generates nine reports designed to identify problem areas and trends in engine 
logistic support, time between engine failures, engine availability, engine readiness degraders, 
and other engine performance factors. The data base used in the ERAP program is the Naval 
Aviation Logistics Data Analysis (NALDA) data base. These reports provide summary 
measures such as mean flight hours between engine failure and other averages and trend lines 
on these averages. They cannot be used, for example, to compute residual reliabilities of 
specified additional lifetimes for engines with known number of repairs and operating times 
since last failure. Nor can they be used to determine optimal maintenance/replacement policies 
for critical engine components. A more detailed analysis of the engine failure data in the 
NALDA data base 1s required to make these types of assessments. Michael R. Caudill [Ref. 1] 
and John A. Malsbury [Ref. 2] made more detailed analyses of TF-34 engine performance in 
their Master of Science theses than those provided under the ERAP program. Mark E. Mlikan 
[Ref. 3] developed an integral equation for determining optimal no-build times for critical 
components of TF-34 engines relative to their hard inspection times. All three of these theses 
used a Weibull distnbution to model the probability distribution of the flight hours between 


engine failures for the TF-34 aircraft engines archived in the NALDA data base. 


Ds NALDA Data Base 

NALDA 1s an operational automated information system. It is a centralized Naval 
Aviation 3M data base supported by the Defense Megacenter, Mechanicsburg, Pa. A user of 
the NALDA data base must make coded queries to the system to obtain specific information. 
This requires considerable experience with NALDA. NALDA contains numerous data 
elements on Navy aircraft engines by engine type and senal number. Anytime an aircraft 
engine 1s inspected, maintained, removed from an aircraft, repaired, transported, reinstalled or 
discarded, a report is submitted that is translated into the NALDA data base. The NALDA 
data base must be “cleansed” before it can be used to perform some detailed statistical 
analyses such as fitting a probability distribution to the empirical cumulative probability 
distribution of operating times between engine failures. M. Caudill did this cleansing of the 
NALDA data for TF-34 engines provided to him by support staff at NAWCADPPE. Details 
of this cleansing process and its purpose are described in Caudill’s thesis [Ref. 1]. 

3. Operating Times Between Engine Failures 

The NALDA data base reports the operating time between engine removals. Engines 
can be removed for non-failure reasons. Caudill [Ref. 1] established a data base of operating 
time between engine failures for TF-34 engines by engine serial number. These times are 
“not” statistically independent because of the manner in which aircraft engines are repaired. 
More often than not, they are not totally overhauled before being reinstalled. Only minor repair 
is made on some occasions. The operating time between the first and second failure of an 
engine is dependent on the operating time to first failure and the extent of the repair. Engines 
with long first failure times and minor first repair will contain many used components with large 
accumulated operating hours when reinstalled. They will likely have shorter operating times to 
the next failure than engines with small operating times to first failure. Consequently, the 


operating times between engine failures, for the same engine, are not statistically independent. 


B. FITTING THE PROBABILITY DISTRIBUTION OF OPERATING TIME 
BETWEEN FAILURES FOR TF-34 ENGINES 


ie Weibull Distribution 

Aircraft engines wear out. Consequently, any probability distribution used to model 
their time to failure should reflect wear-out characteristics. The Weibull distribution has these 
properties. It was used by Caudill [Ref. 1], Malsbury [Ref 2], and Mlikan [Ref. 3] to model 
the distribution of operating times between failures for TF-34 engines. They used a modified 
method of moments procedure to compute the parameters of the fitted Weibull distribution. 
This method used special properties of the Weibull distribution and the exponential distribution; 
the sample mean and the sample standard deviation were used to estimate the parameters of the 
fitted Weibull distribution. It was chosen because it was known to be less sensitive to large 
positive or negative correlations (dependency) in the data than the more standard methods of 
maximum likelihood estimation when statistical independence is assumed. 

2: Research Question 

The analysis performed in this thesis provides a comparison of the accuracies of the 
maximum likelihood estimation procedure (assuming independence) and the modified method 
of moments procedure used by Caudill for fitting a Weibull distnbution to the probability 
distnbution of operating times between engine failures. These comparisons are made for times 
that are randomly generated in a manner that resembles the way in which engine failure times 
occur in the real world. These randomly generated operating times between failures for the 
same engine will not be statistically independent. Once this comparison was completed, it was 
extended to real data using TF-34 engine operating times between failures from the NALDA 
data base. [Ref. 1] 





IU. GENERATION OF FAILURE TIMES 


A. PROCEDURE 


To generate the failure times, a program in S-Plus (See Appendix A) was developed 


with the following properties. 


|. Each engine has 3 components (subsystems). 


2. The failure time distributions of each of these components are Weibull distributions 
with fixed parameters determined a-prion and based on experience. 


3. The data generated is a Type II Censored Data Set, and a total of 3 failure times 
are gathered for each engine. 


With these criteria, the program generates the failure times in the following way: 


Step 1: Three Weibull random numbers are generated. They represent the first failure 
time that each one of the three components could have. These values constitute the 
vector faults]. Three additional Weibull random numbers are generated. They 
represent the second failure time that each one of the three components could have. 
These values constitute the vector faults2.Three additional Weibull random numbers 
are generated. They represent the third failure time that each one of the three 
components could have. These values constitute the vector fau/ts3. 


Step 2: The minimum value of the vector fau/ts/ is selected, this value is the time of 
the first failure of the engine. The program creates a new vector new/fail2 that is the 
difference between the vector fau/ts] and the minimum value selected; it is clear that 
this new vector has a component that is 0. This component is replaced in the vector 
newfail2 by the element of the vector fau/ts2 that has the same index as the first failure. 


Step 3: At this point the engine has two components that have accumulated working 
hours and one that has been just repaired. The vector new/fail2 has two elements that 
have the remainder of the theoretical first failure for these components and one element 
that is the theoretical second failure for the other component. So, again, the minimum 
element of this vector is selected and it is the second failure time of the engine. 


Step 4: The program creates a new vector called newfail3 that is the difference 
between the vector newfai/2 and this mintmum value selected; it is clear that this vector 
has a component that is 0. 


If the second engine failure time has taken place in the same component that had the first 
failure, the 0 element is replaced by the element of the vector fau/ts3 whose index matches that 
of the component that has had the two first failures or by the element of the vector faults2 
whose index matches with the component that had the second failure. Again, there is an 
element in the vector fau/ts3 that is a theoretical failure time for one component just repaired, 
and two elements that have the remainder of theoretical failures. So, once again the minimum 
element of this vector will be the value of the third engine failure time. The vector enginefarls 


gives the three random failure times of the engine. 


B. GENERATION 


An S-Plus program was developed to generate, in the bank of engines, the three failure 
times for the engine. The program can be run for different numbers of engines. In this thesis the 
program was run for 100 engines with a total of 300 failure times and for 20 engines with 60 
failure times. This was done to determine if the number of engines has an apparent effect in the 
results of the study. For each one of these engine sets, 30 replications were generated. These 
replications were exported from S-Plus to Excel with the export tool via the S-Plus file-Export 
menu and used to assess averages on the accuracy of the two methods of estimation - 

The first replication for both 20 and 100 engines can be seen in the Appendix A, all of 


the other replications have been retained on diskettes attached to the master copy of this thesis. 


I. ESTIMATION OF PARAMETERS WITH THE MAXIMUM LIKELIHOOD 
AND THE MODIFIED MOMENT METHOD 
A. PROCEDURE 
The two parameter Weibull distribution analyzed in this thesis has a probability density 
function given by: 


(Oot sco (1) 


or equivalently 


B 
(A ent co (2) 


l 
where 9 = - Customarily B is called the shape parameter and 2 (or 6 ) is called the scale 


parameter. 

To estimate the parameters, the capabilities of the spread sheet program Excel were 
used. The failure times are generated in S-Plus and exported to Excel, where they are displayed 
automatically on the Sheet! Worksheet of Excel. 


The Maximum Likelihood estimators are calculated based on well known formulas 


[Ref. 4]. Solve numerically for B 


, ini: 
KH : E 
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The parameters for the Modified Method of Moments were calculated based on the 
fact that if 7 has a Weibull distribution then 7” has a Exponential distribution with a scale 
parameter of X”. The following demonstration is from [Ref 5]: 


Considering the Weibull distribution given by: 








f(t) = pabihte Oo (6) 
Let x=? so t=x'?, 
dt | 7 
we 
7 aS (7) 
Then 
14 
f(x) =| fF CQ) = o ee) 
l B-l 
lB B2By B o-¥s 
=—x" Br’x” e | (8) 


Therefore: 
X ~ exp(/’). (9) 


It is well known that the mean and standard deviation of the exponential distribution 


are equal. Therefore their statistical estimates « and S should be nearly equal. That is 
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1 < yz _1< 
so] at | DT =e (10) 
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~ fi 


where the 7° are the times between engine failures. 


The solution for 8 in equation (10) will be denoted by 8 . 


l B l n 
Since E(7,°) =6° = (+) , the sample mean, -—»> Z” estimates @ ° , therefore: 
es 


=|" my, (11) 


This is the same estimate obtained for the MLE of 8. 


B. ESTIMATION 


To carry out the estimation of parameters with the two methods, two Macros in Visual 
Basic were developed. One Macro for 20 engines (See Appendix B), that uses 60 failure times, 
and the other for 100 engines (See Appendix C), that uses 300 failure times. 

Basically, these two Macros are similar and therefore they, and procedures for using 
them, can be similarly described. Once the data is exported to Excel, it appears on Sheet]. At 
this point it is necessary to open the file with the Macro, called Macro20.x/s for the sample of 
20 engines and Macro/00.xls for the sample of 100 engines respectively. Next the window 


with the Macro is hidden from the same Sheet/.Then the Macro can be run. To run the Macro, 


choose Macro from the Tools menu. From Macro/Name Reference choose the name of the 
Macro and click on run. The Macro does all the previous calculations and operations needed 
to estimate the parameters and to entitle the entire Worksheet. After these calculations the 
Macro stops. At this point it 1s necessary to introduce the following data: 

In Cell B6 introduce the value previously calculated by the macro for the initial 
value of the shape parameter which ts in Cell B4. 

This initial parameter was introduced to assure the convergence with the Tool Goal 
Seeker, because without it, a solution may not be reached. The initial value is given by [Ref. 


6]: 
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Using these choices for the initial shape parameter, an average of only 3.5 Newton 


iterations were needed to achieve four-place accuracy in the solution for 8. 

After introducing the initial value given in Cell B4 into Cell B6, cell C6 must be 
highlighted. Then the Tools Menu is selected, Goal Seek 1s selected, a 0 1s placed in the “To 
Value,” and B6 is entered as the “Changing Cell.” 


Equation (3) 1s in Cell C6 and the Goal Seeker searches for the value of § which 
satisfies that equation. There were no non-convergence problems in any of the 70 
replications made with that formula. Once the solution 1s achieved, the final value for the 
shape parameter appears in Cell B6. The solutions for 24 and 9 appear in Cell B9 
(Equation (4)) and B11 (Equation (5)) respectively. These solutions for #, 2 and 0 are 
the maximum likelihood estimates for these parameters of the Weibull distribution that 1s 
used to model the probability distribution of the family of times between failures on all 
engines. These values are printed on the left side of Sheet2. 

On the right side of Sheet 2 are the calculations carried out by the Modified 


Moment Method. 


There were convergence problems in calculating 8 using Equation (10). Amazingly 
the same formula may converge or not depending on the terms on each side. Through tral and 
error it was discovered that by wniting the formula in the following way it was possible to get 
convergence in all replications: 

eg 
n = i=l 
sal Pe 0.5. (13) 
pa (Ti? wy | 


f=) 


In Cell Q6 it is necessary to introduce the initial shape parameter given in Cell Q4 
(Equation (12)). Once done, highlight Cell P7 (Equation (13)), select the Tools Menu, select 
Goal Seek and depending on the number of engines, enter in “To Value” 7.811334 (20 
engines, 60 failure times, 60/ 59 ) or 17.349447 (100 engines, 300 failure times, 300/ 299 ). 
In the “Changing Cell, Q6,” the Goal Seeker looks for the B value that will satisfy Equation 
(13). Cells Q9 and Q11 contain the values for the scale parameter 4 and 8 determined by 
Equation (11) and (5). 

A total of 30 replications for each set of engines were done. Afterwards, for the 
reasons explained later, 10 more replications were done using real data. 

Printouts of the failure times and parameter calculations for the 20 and 100 engines first 
replication are provided in Appendix E. 


A summary of the parameters obtained is shown in the following tables: 


1] 


Table 1. Twenty Engines = 60 Failure Times 


20 Engines = 60 Failure Times 


Maximum Likelihood Estimator 
Shape Scale 







41 0.973000000 | _251.70530000 
3] 1.105104390 | 279.48190000 

242.34360000 

282.57760000 

272.41390000 

277 44380000 
F401 1.108272000 | 263.760200000) 
2| 1.358459000 |_ 245.95380000 
28103880000 
4| 1.054978965 |_223.27680000 
6| 1.228101693 | _264.36450000 
243.73460000 
47| 1.355617000 | _268.067900000! 
48| 1.308393000 | 309.037200000! 
49 
'4| 1.336786000 |_ 271.132800000 
‘231 1.200747000 | _288.213200000 
'25| 1.400818000 |_310.907900000! 
'26| 1.231510000 |_277.128900000) 
'29| 1.295133208 | 246.789500000 
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Modified Moment Method 
Shape Scale 








273.090800000 


O jO 





Average 1.203679065 272.060400000 
Std.Dev. 0.145886716  25.376732459 


1.235443638 274.848080000 
0.133749260 25.937151577 
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Replication Number 


Average 
Std. Dev. 


oO 


“21 4.147177000| 258.38230000 
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Table 2. One-Hundred Engines = 300 Failure Times 


100 Engines = 300 Failure Times 


Maximum Likelihood Estimator 
Shape 


Scale 


1.261028000{ 292.55690000 
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1.19224900 
1.19139300 
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10016700 
13262100 
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26850900 
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274659000! 270.850700000 
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Modified Moment Method 
Shape Scale 


1.288232000| 294 659500000 
1.182793000] 261.24640000 
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1271795400 278.553426667 
0.070516179 13.009908323 
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IV. ACCURACY ANALYSIS 


A. PROCEDURE 


The accuracy of each method was measured using the Kolmogorov-Smimov goodness 
of fit Test and the S-Plus software package. The null hypothesis in S-Plus is “the True cdf 
(empirical) equals the Weibull distribution”. This test has the advantage that it uses continuous 
populations. Therefore it is not necessary to group the data. 

This test has the disadvantage that when distribution parameters are estimated from the 
same sample of data, the test produces a conservative result; that is, the probability of a Type I 
error will be smaller than specified [Ref 6]. This drawback has only a relative importance in the 
present study. The K-S Test p-value is a conservative result whether the parameters are 
estimated by the Maximum Likelihood Method or the Modified Moment Method. 

A S-Plus function called Zocompare (See Appendix E) was created to calculate the K- 
S statistics for the set of parameters calculated in each replication using both methods. This S- 
Plus function computes the K-S Test Statistic (Supremum) and the p-value for each method, 
and it also plots the empirical and hypothesized distributions. In Appendix E the reader will find 


the output of the function for the first replication of each set of parameters. 


B. ACCURACY 


The calculations for all replications and the values obtained for each set of engines are 
presented in Tables 3 and 4. The letter “ B” marks the method that achieved a better fit, using 
as reference the K-S statistic; that is, the method with the smaller K-S statistic 1s the more 
accurate procedure for that replication. 

For the set of 20 engines, the MLE method gave a better fit in 17 of the 30 replications. 
For the set of 100 engines, the MLE method gave a better fit in 20 of the 30 replications. 


Is 


Replication Number 


shape 


























0.973000000 
41.105104390 
H 1.052707000 
44 1.005354000 
1.211857000 
1.523761000 
4.053815000 
1.015065000 
0.961137000 
1.108272000 
1.258861000 
1.358459000 
1.323305000 
1.054978965 
1.228101693 
1.161814103 
4.35561 7000 
84 1.308393000 
1.391697000 
1.141031467 
4.336786000 
41.389460000 
1.200747000 
1.29742 7000 
1.400818000 
1.231510000 
4.048038281 
1.093433000 
1.295133208 
1.224688831 
mean 
1.203679065 

stdev 
0.145886716 








Table 3. Better Fit Twenty Engines 


20 Engines = 60 Failures times 


M.L.E. 


scale 


251.705300000 
279.481900000 
242.343600000 
282.577600000 
272.413900000 
277.443800000 
243.377800000 
330.527000000 
251.822700000 
263. 760200000 
317.378200000 
245.953800000 
281.038800000 
223.276800000 
264.364500000 
243.734600000 
268.067900000 
309.037200000 
284.440700000 
242.289400000 
271.132800000 
244.898600000 
288.213200000 
293.585500000 
310.907900000 
277.128900000 
277.331100000 
273.733700000 
246.789500000 
303.055100000 
mean 
272.060400000 

stdev 
25.376732459 


k-s 


0.054600000 
0.088300000 
0.076300000 
0.061700000 
0.101800000 
0.086700000 
0.09 1500000 
0.052700000 
0.088000000 
0.064500000 
0.101500000 
0.104300000 
0.075700000 
0.088000000 
0.070300000 
0.078800000 
0.070200000 
0.049400000 
0.083700000 
0.091800000 
0.107000000 
0.085000000 
0.083700000 
0.058400000 
0.071100000 
0.062500000 
0.078500000 
0.078000000 
0.055600000 
0.083900000 
mean 
0.078116667 

stdev 
0.015610639 





p-value 


0.994000000 
0.737300000 
0.876000000 
0.976300000 
0.562800000 
0.757300000 
0.696500000 
0.996200000 
0.742000000 
0.964400000 
0.566300000 
0.531000000 
0.881800000 
0.741900000 
0.928000000 
0.850400000 
0.929200000 
0.998600000 
0.794200000 
0.692600000 
0.498400000 
0.778500000 
0.795100000 
0.986700000 
0.921800000 
0.973200000 
0.853300000 
0.859200000 
0.992500000 
0.792200000 
mean 
0.822256667 








stdev 


0.146004782 








shape 





0.985725000 
1.129925000 
1.036553000 
1.079900000 
1.102670000 
1.442500000 
1.309674000 
1.006751000 
1.078819000 
1.084296000 
1.395249000 
1.211490000 
1.461631834 
1.249979811 
1.182801000 
1.239576000 
1.36 7488000 
1.35 7252000 
1.406588000 
1.32 4936416 
1.223524000 
1.280865198 
1.263844000 
1.260780021 
1.325062000 
1.320218000 
1.043436000 
1.248177000 
1.240951 783 
1.402646081 
mean 
1.235443638 







































0.133749260 
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stdev 


M.M.M. 


scale 


253.139700000 
281.796500000 
240.792500000 
291.030800000 
263.279100000 
273.090800000 
262 .946000000 
329 .425200000 
264.611000000 
261.569300000 
328.224900000 
237.176300000 
289.757400000 
236.965600000 
260.845800000 
249.537300000 
268.798900000 
312 .609300000 
285 .359900000 
255.330700000 
263.236700000 
238.618200000 
293.349300000 
290.891900000 
305 .455500000 
283.439100000 
276.830600000 
286 .965800000 
243 .300600000 
317.067 700000 
mean 








k-S 


0.050900000 
0.08 1000000 
0.073000000 
0.086100000 
0.097300000 
0.08 1000000 
0.115700000 
0.053100000 
0.098600000 
0.06 1500000 
0.126600000 
0.121500000 
0.110800000 
0.087200000 
0.065700000 
0.101500000 
0.072500000 
0.054400000 
0.083400000 
0.116800000 
0.07 7500000 
0.092 700000 
0.096600000 
0.059800000 
0.071100000 
0.057400000 
0.077100000 
0.084100000 
0.06 7600000 
0.095900000 
mean 








p-value 








0.997700000 
0.826600000 
0.906600000 
0.765600000 
0.62 1600000 
0.826400000 
0.398400000 
0.995900000 
0.604700000 
0.977100000 
0.290900000 
0.338000000 
0.452600000 
0.751700000 
0.958200000 
0.567200000 
0.910600000 
0.994300000 
0.798500000 
0.386000000 
0.863800000 
0.681600000 
0.629800000 
0.982700000 
0.922100000 
0.989000000 
0.86 7800000 
0.790390000 
0.947000000 
0.639300000 
mean 


274.848080000 0.083946667 0.756066667 


stdev 


stdev 


25.937151577 0.020788358 0.214131410 


stdev 


Replication Number 





Omwsh th R Oh = 


shape 


4.261028000 
1.147177000 
1.192249000 
4.191393000 
4.100167000 
41.132621000 
4.220806000 
4.206135000 
4.254734000 
1.284450000 
1.277785000 
1.274025000 
4.192249000 
1.194163000 
1.216995000 
1.169718000 
1.283386000 
1.416971000 
1.264662000 
1.342221000 
1.229806000 
4.199420000 
1.189326000 
4.199541000 
1.268509000 
1.257310000 
1.274659000 
4.208075000 
1.220756000 
4.127942000 
mean 
1.226599300 
standev 

0 063780678 


Table 4. Better Fit One-Hundred Engines 


100 Engines = 300 Failures times 


M.L.E. 


scale 


292.556900000 
258.382300000 
267.935 700000 
307.254500000 
265.867 700000 
273.303500000 
282 .032200000 
267.01 7000000 
291.692500000 
285.463200000 
285.850800000 
263.121900000 
267.935 700000 
271.996400000 
272.494100000 
293.370600000 
255.640200000 
294.799600000 
268.682800000 
279.712400000 
287.109500000 
267.635600000 
288.996500000 
263.550100000 
261.778100000 
256.670600000 
270.850700000 
276.652800000 
264.013900000 
279.203900000 
mean 
275.185723333 
standev 
12.379509005 


k-s 


0.03 1900000 
0.027400000 
0.047000000 
0.055800000 
0.060600000 
0.043200000 
0.034300000 
0.037900000 
0.02 71600000 
0.022 700000 
0.044600000 
0.047400000 
0.029400000 
0.039600000 
0.03 1400000 
0.034800000 
0.03 1800000 
0.039800000 
0.053800000 
0.039700000 
0.041300000 
0.036000000 
0.032800000 
0.032000000 
0.028200000 
0.026400000 
0.038700000 
0.035200000 
0.033800000 
0.035 700000 
mean 
0.037760000 
standev 
0.009086679 


p-value 


0.920600000 
0.977600000 
0.522300000 
0.306900000 
0.220000000 
0.630900000 
0.873100000 
0.782600000 
0.999000000 
0.997900000 
0.589000000 
0.511700000 
0.957700000 
0.735300000 
0.928900000 
0.861700000 
0.922000000 
0.729300000 
0.349300000 
0.730800000 
0.685900000 
0.830800000 
0.903200000 
0.919200000 
0.971000000 
0.984800000 
0.759500000 
0.851300000 
0.882800000 
0.838800000 
mean 
0.772463333 
standev 
0.209187169 
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shape 


1.288232000 
1.182793000 
1.250620000 
1.243339000 
1.187671000 
1.205131000 
1.334481000 
1.246956000 
1.271828000 
1.292318000 
4.387920000 
1.341512000 
1.250620000 
1.228539000 
1.231497000 
1.200586000 
1.259945000 
1.422087000 
1.379585000 
1.448534000 
1.26996 7000 
41.281756000 
1.207016000 
1.346121000 
1.201983000 
1.249969000 
4.311836000 
4.237968000 
1.219361000 
1.173691000 
mean 
1.271795400 
standev 
0.070516179 


M.M.M. 


scale 


294.659500000 
261.246400000 
272.423600000 
305.757700000 
272.564100000 
279.422100000 
290.504100000 
270.084700000 
293.041800000 
286.039200000 
293.508300000 
267.591 700000 
272.423600000 
274.677200000 
273.612100000 
296.122500000 
254.068200000 
295.112700000 
276.306900000 
286.372700000 
290.254800000 
273.804000000 
290.494900000 
273.747200000 
256 .990000000 
256.167900000 
273.423600000 
278.973700000 
263.912300000 
283.295300000 
mean 
278.553426667 
standev 
13.009908323 


0.038200000 
0.032 700000 
0.051000000 
0.046600000 
0.045600000 
0.042500000 
0.048900000 
0.036300000 
0.025500000 
0.020800000 
0.058300000 
0.048900000 
0.035300000 
0.048600000 
0.029400000 
0.041300000 
0.033000000 
0.038700000 
0.065700000 
0.06 7500000 
0.040300000 
0.044400000 
0.033600000 
0.048900000 
0.041500000 
0.027 700000 
0.048400000 
0.028200000 
0.034200000 
0.031400000 
mean 
0.044113333 
standev 
0.010947594 


p-value 


0.774000000 
0.905100000 
0.416900000 
0.532800000 
0.567900000 
0.6517400000 
0.469300000 
0.823000000 
0.989900000 
0.999500000 
0.260200000 
0.469900000 
0.848500000 
0.477100000 
0.957500000 
0.685300000 
0.898900000 
0.759500000 
0.150500000 
0.129900000 
0.713200000 
0.596300000 
0.886800000 
0.469200000 
0.680900000 
0.975800000 
0.482200000 
0.971300000 
0.874100000 
0.928300000 
mean 
0.677973333 
standev 
0.247541499 


If the MLE and MMM methods are equally good then one would expect that the 
probability the MLE gives a better fit ( using the K-S statistic) should be 0.5. Let 2 represent 
the probability the MLE gives a better fit; the above sets of data can then be used to test Ho: 2 
= 0.5 versus H,: m #0.5. For the set of 20 engines, 17 successes were observed out of 30 
replications ( MLE gave a smaller K-S value than MMM); the p-value for testing Hp : 2 = 0.5 
is 0.5847 so Hp would be accepted for any « < 0.5847. For the set of 100 engines, 20 of the 
30 replications gave successes; in this case the p-value 1s 0.0987 and Hp is rejected for any a 
> 0.0987. Certainly this comparison of MLE and MMM, using the binomial test, shows no real 
difference for 20 engines and marginally indicates supenionty of MLE for 100 engines. 

If the K-S statistic mean and standard deviation are analyzed it is clear that the MLE. 
gives the smaller mean and standard deviation of the two cases. 

For this data, it seems that despite the fact that there is some dependence between the 
failure times, the MMM does not give a better fit than the MLE. This is true at least for the 
failure times generated with the simulation program used in this study. This is likely due to the 
fact that the failure times for different engines (among engines) are independent. 

At this point real failure data was examined. This real data is from [Ref: 1], and from 
the cleansed data collected in Appendix B (Five-Year Time Between Failures Data Base). This 
data was cleansed, removing zeros as failure times. Ten samples of data were taken from the 
data set of operating times between failures on TF 34 aircraft engines. Five samples were of 60 
failures times and five samples were of 300 failure times. 

The same procedure to estimate the parameters and analyze the accuracy of the results 
were applied. The results are shown in Tables 5 and 6. 

With the real data the result is the same, a better fit is obtained more often with the 
MLE method than with the MMM method. In the 60 failure times case the MLE method was 
the better procedure in four of the five samples, and in the 300 failure times case it was better 


three of the five samples using the K-S statistic as measure of accuracy. 
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It is important to point out that in these real data, the engine structure is unknown. 
There is no information about how many elements the engine has, or of course the way that 


these components failed. 


Table 5. Estimate Parameters and Analyze Accuracy Twenty Engines 


20 Engines = 60 Failures times 


M.L.E. M.M.M. 
shape scale k-s p-value shape scale K-s p-value 









-441.554927 651.243400 0.103000 0.547400 B } 1.732451 669.495500 0.103200 0.544600 
211.435165 710.748800 0.079700 0.840000 B [1.735014 748.435600 0.130200 0.261300 
3] 1.531311 603.432400 0.072600 0.909800 1.635563 613.677100 0.051800 0.997100 B 
41 1.436419 609.475100 0.103800 0.537100 B {1.485100 615.440400 0.114900 0.406200 
$1 1.535733 579.331200 0.135400 0.221300 B {1.490634 574.988000 0.145900 0.155400 
mean mean mean mean mean mean mean mean 
1.414134 642.431620 0.048360 0.511380 411.587754 663.677300 0.057660 0.289300 1 


Replication 
Number 


stdev stdev stdev stdev stdev stdev stdev stdev 
0.045711 14.890226 0.008406 0.185231 0.023623 16.834529 0.006280 0.138620 


C. ANALYSIS OF THE RESULTS 


The author expected “ a-prion” a better fit with the MMM than with the MLE, based 
on the fact that the generated operating times between failures for the same engine are not 
statistically independent. But the resulting fit for the data was better with the MLE more often 
than with the MMM. This result is further analyzed here, by examining how the dependence 
between failure times affect each method. 

To assess dependence between failure times, the correlation and the grade of 
dependence between the failure times was calculated for each replication. The correlation 
coefficients were calculated and added at the end of each replication (See Appendix A). In 
addition the Kendall-Tau independence test was performed. The null hypothesis in S-Plus for 
this test is “ Tau is equal to 0” that is the two samples are independent. This analysis was 
performed for all possible pairwise combinations of the failure times of an engine; namely first- 
second, first-third, second-third correlations, tau coefficients and p-values were calculated. For 


the Kendall-Tau test the minimum, maximum and mean p-values of those pairwise 
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combinations of each replication were added to the tables. The results are shown in Tables 7 
through 10. 


Table 6. Estimate Parameters and Analyze Accuracy One-Hundred Engines 


100 Engines = 300 Failures times 
M.L.E. M.M.M. 
shape scale k-s p-value shape scale k-s p-value 














1.579829 632.033000 0.045700 0.559000 B 
1.586168 681.405000 0.063000 0.184300 
1.564422 666.063600 0.057300 0.279000 
1.575530 665.621100 0.061900 0.200100 B 
1.632821 673.263800 0.060400 0.224100 
mean mean mean mean 
1.587754 663.677300 0.057660 0.289300 2 
stdev stdev Stdev stdev 
0.023623 16.834529 0.006280 0.138620 


44 1.462681 618.916600 0.047200 0.515100 
“74 1.439461 663.159800 0.045200 0.571700 B 
~341.366409 641.062700 0.040500 0.709200 B 
~A11.351920 636.911100 0.064600 0.163100 

541.450201 652.107900 0.044300 0.597800 
mean mean mean mean 
1.414134 642.431620 0.048360 0.511380 3 

stdev stdev stdev stdev 
0.045711 14.890226 0.008406 0.185231 
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From those tables may be deduced that with the set of 20 engines the correlation 1s 
clearly larger than with 100 engines. When the independence is tested with the Kendall-Tau 
test, the values of the tau’s coefficients are similar to the correlation coefficients, and have 
larger values for the set of 20 engines than the set of 100 engines. In general for this set of data, 
the correlation coefficients are small and the p-values in the Kendall-Tau are not significant 
enough to reject the null hypothesis of independence at the usual values. That means the 
dependence obtained in the simulation is small, specially with the set of 100 engines; this may 
be unfavorable to the performance of the MMM. 

To summanize the results, Figures 1 through 4 present the min and max Kendall-Tau 
test p-value for each replication for the method that gives the best fit. 

From these figures and specially from the max p-values, it may be deduced that, for this 
data, in general the larger the max p-value (more independence) the better the fit with the MLE 
method. Conversely the smaller the max p-value (less independence) the better the fit with the 


MMM method. 
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Table 7. Correlation - 20 Engines 
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Table 8. Correlation - 100 Engines 
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Table 9. Kendall -Tau Test - 20 Engines 
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Table 10. Kendall -Tau Test - 100 Engines 
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Figure 1. Maximum p-Value for 20 Engines 
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Figure 2. Minimum p-Value for 20 Engines 
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Figure 3. Maximum p-Value for 100 Engines 
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Figure 4. Minimum p-Value for 100 Engines 
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V. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY 


The analytical approach taken to develop this thesis and topics covered has been the 
following: 
1. Develop a program in S-Plus to generate the failure times for an engine with 3 


subsystems, where the distribution of the first failure times of these subsystems is 
Weibull and the engine 1s not restored to new condition after each engine failure. 


2. Estimate the parameters for the Weibull distribution that fits the generated engine 
failure times with the Maximum Likelihood Estimator and a Modified Method of 
the Moments using commercial spread-sheet software. 


3. Measure the accuracy of these two methods with the Kolmogorov-Smimov Test. 
4. Apply the same analysis to a given set of real engine failure times. 


5. Compare the accuracy of the two methods, and discuss reasons for their 
differences. 


B. CONCLUSIONS 


The two methods in general give good results when they are applied to fit a Weibull 
distribution to the failure times generated in the simulation. The resulting fit for these failure 
times was better with the Maximum Likelihood Estimator more often than with the Modified 
Moment Method. For this reason, and the fact that the MLE method has superior theoretical 
properties in the presence of independence, the author recommends the MLE method to 
estimate the parameters for scenarios similar to those simulated in this Thesis. 


With these conclusions it is important to point out the following considerations: 


2g 


The MMM has the advantage that fewer calculations are necessary to get the results. 
In the case of the Macros developed in this study, only the initial values and two columns of 
calculations are necessary. Four columns of calculations are necessary in the MLE case. 

There was a wide range of values for the parameters of the fitted distributions across all 
the replications. This is displayed in Table 11. 


Table 11. Parameters Range 


20 Engines 
MLE MMM 
0.9611 <= shape <= 1.5238 0.9857 <= shape <= 1.4616 
223.28 == Scale —— 353 236.97 <= scale <= 329.43 


100 Engines 
MLE MMM 
1.1002 <= shape <= 1,3422 1.1737 <= shape <= 1.4485 
255.64 <= scale <= 301.25 254.07 <= scale <= 305.76 





The dependence obtained in the simulation is not as significant as was expected. 


Ce RECOMMENDATIONS 


These two methods should be studied when greater dependence exists among the data 


and when the dependency is introduced in a different manner than used in this thesis. 


APPENDIX A. PROGRAM TO GENERATE FAILURE TIMES 


function(n) 
t 
# generate 
# This function generates times randomly in a manner that resembles the way in which 
# engines failures occur in the real world. 

totalreplic <- NULL 

for(count in 1:n) { 
# Step | 
# Three Weibull random numbers are generated, each one represents the first failure time 
# that each one of the components could have. These values constitute the vector faults]. 
# Three Weibull random numbers are generated, each one represents the second failure time 
# that each one of the components could have. These values constitute the vector faults2. 
# Three Weibull random numbers are generated, each one represents the third failure time 
# that each one of the components could have. These values constitute the vector faults3. 

x11 <- rweibull(1, 1.2, 900) 

x12 <- rweibull(1, 1.3, 800) 

x13 <- rweibull(1, 1.4, 600) 

faults] <- c(x]l1, x12, x13) 

x21 <- rweibull(1, 1.2, 900) 

x22 <- rweibull(1, 1.3, 800) 

x23 <- rweibull(1, 1.4, 600) 

faults2 <- c(x21, x22, x23) 

x31 <- rweibull(1, 1.2, 900) 

x32 <- rweibull(1, 1.3, 800) 

x33 <- rweibull(1, 1.4, 600) 

faults3 <- c(x31, x32, x33) 

failures <- matrix(c(faults1, faults2, faults3), ncol = 3) 
# Step 2 
# The minimum value of the vector faults1 1s selected, this value is the time of the first 
# failure of the engine. The program creates a new vector newfail2 that is the difference 
# between the vector faults] and the minimum value selected, it is clear that this new vector 
# has a component that is 0, this component is replaced in the vector newfail2 for the 
# element of the vector faults2 that match with the component which had the first failure. 
# Choose the first failure time, that is the minimum 
# time from the vector faults] 

fl <- min(faults1) 
# the next step is to choose the second failure time, to make this first 
# the above value min(faults1) must be substituted for a new time 
# generated in the vector faults2 and precisely in the same component. 
# After this the new minimum is selected 

newfail2 <- (faults! - f1) 

11 <- match(0. newfail2) 

newfail2 [11] <- faults2 [1] 


# Step 3 


# At this point the engine has two components that have accumulated working hours and 
# one that has been just repaired. The vector newfail2 has two elements that have the 


3] 


# remainder of the theoretical first failure for these components and one element that is the 
# theoretical second failure for the other component. So, again, the minimum element of 
# this vector is selected and it is the second failure time of the engine. 


f2 <- min(newfail2) 


# Step 4 

# The program creates a new vector called newfail3 that is the difference between the vector 
# newfail2 and this minimum value selected, it is clear that this vector has a component that 
# 1s 0. 


# Choose the third failure time, here there are two possible 
# cases, that the second failure is in the same component that 
# the first one or not, these two condition are programmed in this order 
# 
newfail3 <- (newfail2 - f2) 
ili <- match(O, newfail3) 
if(ii = ill) £ 
newfail3 [11] <- faults3 [11] 
f3 <- min(newfail3) 


} 
else if(ii '!= iii) { 
newfail3 [iii] <- faults2([111] 
f3 <- min(newfail3) 
} 
enginefails <- c(fl, f2, f3) 
#cumenginefails <- cumsum(enginefails) 
#engfails <- matnx(c(failures, enginefails, cumenginefails), byrow = 
# T,ncol = 3) 
totalreplic <- rbind(totalreplic, (enginefails)) 
} 
totalreplic 
} 
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Row Names V1 V2 V3 
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APPENDIX B. MACRO 20 ENGINES 


"MACRO FOR 20 ENGINES 
" thesis Macro 
‘Macro recorded by FERNANDO GUIMIL 


Sub thesis() 
Sheets. Add 
Sheets("Sheet2").Select 
Range("A1:C2").Select 
ActiveCell.FormulaR1C1 = "MAXIMUM LIKELIHOOD ESTIMATOR" 
Range("A1:C2").Select 
Range("A2"). Activate 
Selection.Font.Bold = Tre 
With Selection 
Horizontal Alignment = x]Center 
.VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
With Selection 
Horizontal Alignment = xlCenterAcrossSelection 
Vertical Alignment = xlBottom 
.WrapText = True 
Orientation = xlHorizontal 
End With 
Range("D1").Select 
ActiveWindow.SmallScroll ToRight:=10 
With Toolbars(8) 
Left = 241 
.Top = 228 
End With 
ActiveWindow.LargeScroll ToRight:=-1 
Range("P1:R1").Select 
ActiveCell.FormulaR1C1 = "MODIFIED MOMENT METHOD" 
Range("P1:R1").Select 
Range("Q1"). Activate 
Selection.Font.Bold = True 
With Selection 
.HorizontalAlignment = xiCenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
With Selection 
-Honzontal Alignment = xIlCenterAcrossSelection 
. VerticalAlignment = xlBottom 
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.WrapText = True 
Orientation = xlHorizontal 
End With 
Range("N1").Select 
ActiveWindow.LargeScroll ToRight:=-] 
With Toolbars(8) 
Left = 436 
Top = 226 
End With 
Range("A2").Select 
ActiveCell.FormulaR1C1 = "shape parameter" 
With ActiveCell. Characters(Start:=1, Length:=15).Font 
.Name = "Anal" 
FontStyle = "Bold" 
size = 10 
Strikethrough = False 
Superscript = False 
Subscript = False 
OutlineFont = False 
Shadow = False 
_Underline = xINone 
ColorIndex = xlAutomatic 
End With 
Range("B3").Select 
ActiveCell.FormulaR1C1 = "initial beta" 
Range("B3").Select 
Selection.Font.Bold = True 
With Selection 
-HorizontalAlignment = xlCenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
Range("B5"). Select 
ActiveCell.FormulaR1C1 = "final beta" 
Range("B5").Select 
Selection.Font.Bold = True 
With Selection 
-HorizontalAlignment = x!Center 
.VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
Range("A7:C7"). Select 
ActiveCell.FormulaR1IC1 = "scale parameter" 
Range("A7:C7"). Select 
Range("B7"). Activate 
Selection. Font.Bold = True 
With Selection 
-HorizontalAlignment = xICenter 
.VerticalAlignment = xIBottom 
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.WrapText = False 

Orientation = xiHorizontal 
End With 
Range("A7").Select 
ActiveCell.FormulaR1C1] = "" 
Range("B7").Select 
ActiveCell.FormulaR1C1 = "scale parameter" 
Range("B8").Select 
ActiveCell.FormulaR1C1] = “landa" 
Range("B8").Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xICenter 

.VerticalAlignment = xlBottom 

.WrapText = False 

Orientation = xlHorizontal 
End With 
ActiveWindow.SmallScroll ToRight:=0 
Range("A2:C2").Select 
Selection.Copy 
ActiveWindow.SmallScroll ToRight:=] } 
Range("Q2").Select 
ActiveSheet.Paste 
Application.CutCopyMode = False 
Selection.Clear 
ActiveWindow.LargeScroll ToRight:=-2 
Range("A2:C2").Select 
Selection.Copy 
Active Window.SmallScroll ToRight:=9 
Range("P2").Select 
ActiveSheet.Paste 
Range("O4").Select 
Active Window.LargeScroll ToRight:=-1 
Range("B3"). Select 
Application.CutCopyMode = False 
Selection.Copy 
ActiveWindow.SmaliScroll ToRight:=11 
Range("Q3").Select 
ActiveSheet.Paste 
ActiveWindow.LargeScroll ToRight:=-2 
Range("BS5S").Select 
Application.CutCopyMode = False 
Selection. Copy 
ActiveWindow.SmallScroll ToRight:=]2 
Range("Q5").Select 
Application.CutCopyMode = False 
Selection.Copy 
ActiveSheet.Paste 
Range("R7").Select 
ActiveWindow.LargeScroll ToRight:=-2 
Range("B5").Select 
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Application. CutCopyMode = False 
Selection.Copy 
ActiveWindow.SmallScroll ToRight:=12 
Range("Q5"). Select 
ActiveSheet.Paste 
ActiveWindow.LargeScroll ToRight:=-2 
Range("A7:C8"). Select 
Application. CutCopyMode = False 
Selection.Copy 
ActiveWindow.SmallScroll ToRight:=1 1 
Range("P7"). Select 
ActiveSheet.Paste 
Range("N7").Select 
ActiveWindow.LargeScroll ToRight:=-2 
Range("D11").Select 
Application. CutCopyMode = False 
Selection.Copy 
Range("D7"). Select 
ActiveSheet.Paste 
Range("B8").Select 
ActiveWindow.SmallScroll ToRight:=13 
ActiveWindow.LargeScroll ToRight:=-2 
ActiveWindow.SmallScroll ToRight:=14 
ActiveWindow.LargeScroll ToRight:=-2 
Range("E1").Select 
Application.CutCopyMode = False 
ActiveCell.FormulaR1Cl = "Number of observations" 
Range("E1").Select 
Selection.Font.Bold = True 
With Selection 
-HorizontalAlignment = x]Center 
.VerticalAlignment = xlBottom 
.WrapText = False 
.Onentation = xlHorizontal 
End With 
With Selection 
-HonzontalAlignment = x!CenterAcrossSelection 
. Vertical Alignment = xIBottom 
.WrapText = True 
Orientation = xlHorizontal 
End With 
Range("E2"). Select 
ActiveCell.FormulaR1C1 = "n" 
Range("E2").Select 
With Toolbars(8) 
Left = 536 
lop = 23 
End With 
Sclection.Font.Bold = True 
With Selection 
-HorizontalAlignment = xlCenter 
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-VerticalAhgnment = xIBottom 

.WrapText = False 

.Onentation = xlHorizontal 
End With 
Sheets("Sheet1").Select 
Range("B2:B21").Select 
Selection. Copy 
Sheets("Sheet2").Select 
Range("D10").Select 
ActiveSheet.Paste 
Sheets("Sheet1").Select 
Range("C2:C21").Select 
Application.CutCopyMode = False 
Selection.Copy 
Sheets(""Sheet2").Select 
Range("E10").Select 
ActiveSheet. Paste 
Sheets("Sheet 1").Select 
Range("D2:D21").Select 
Application.CutCopyMode = False 
Selection.Copy 
Sheets("Sheet2").Select 
Range("F10").Select 
Application.CutCopyMode = False 
Selection.Copy 
Sheets("Sheet 1").Select 
Application.CutCopyMode = False 
Selection.Copy 
Sheets("Sheet2"). Select 
ActiveSheet.Paste 
Range("D10:D29"). Select 
Application.CutCopyMode = False 
Selection.Copy 
Range("G10").Select 
ActiveSheet.Paste 
Range("E10:E29"). Select 
Application. CutCopyMode = False 
Selection.Copy 
ActiveWindow.SmalScroll Down:=16 
Range(""G30").Select 
ActiveSheet.Paste 
ActiveWindow.LargeScroll Down:=-1 
Range("F10:F29"). Select 
Application.CutCopyMode = False 
Selection.Copy 
ActiveWindow.SmalScroll Down:=35 
Range("G50").Select 
ActiveSheet.Paste 
Range("D42").Select 
ActiveWindow.LargeScroll Down:=-2 
Range("G10:G69"). Select 
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Application.CutCopyMode = False 
Selection. Copy 
Range("H10").Select 
ActiveSheet.Paste 
Application.CutCopyMode = False 
Selection.Sort Key! :=Range(""H10"), Order! :=xlAscending, Header:= _ 
xIlGuess, OrderCustom:=1, MatchCase:=False, Orientation:= _ 
xlITopToBottom 
Range("B10").Select 
ActiveCell.FormulaR!C1 = "theta" 
Range("B10").Select 
Selection.Font.Bold = True 
With Selection 
Horizontal Alignment = xICenter 
.VerticalAlignment = xIBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
Range("E3").Select 
ActiveCell.FormulaR1C1 = "60" 
Range("E3").Select 
Selection.Font.Bold = True 
With Selection 
-HorizontalAlignment = x!Center 
.VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
Active Window. ScrollRow = 1 
Range("G1").Select 
ActiveCell.FormulaR!C1 = "Failures times Ti" 
Range("G1").Select 
With Selection 
-HorizontalAlignment = xlCenterAcrossSelection 
.VerticalAlignment = xlBottom 
.WrapText = True 
.Onentation = xlHorizontal 
End With 
Selection.Font.Bold = True 
Range("G9").Select 
ActiveCell.FormulaRiCl = "Ti" 
Range("G9"). Select 
Selection. Font.Bold = True 
With Selection 
-HorizontalAlignment = x]Center 
.VerticalAlignment = xiBottom 
.WrapTent = False 
Orientation = xlHorizontal 
End With 
Selection. Copy 
Range("H9").Select 


ActiveSheet.Paste 
Range("G9").Select 
Application.CutCopyMode = False 
Selection.Clear 
ActiveWindow.SmallScroll ToRight:=3 
Range("19"). Select 
ActiveCell.FormulaR!1C1 = "Ln Ti" 
Range("19").Select 
Selection.Font.Bold = True 
With Selection 
-HorizontalAlignment = xICenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
Range("I10").Select 
ActiveCell.FormulaR1C1 = "=LN(RC[-1])" 
Range("I10").Select 
Selection. AutoFill Destination:=Range("110:169"), Type:= _ 
xIFillDefault 
Range("I10:169"). Select 
Range("H70").Select 
ActiveCell.FormulaR1C1 = "SUM" 
Range("H70").Select 
Selection.Font.Bold = True 
With Selection 
.HorizontalAlignment = xICenter 
. VerticalAlignment = xIBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
Selection.Copy 
Range("170"). Select 
ActiveSheet.Paste 
Range("H71").Select 
Application.CutCopyMode = False 
ActiveCell.FormulaR1C1] = "=SUM(R[-61]C:R[-2]C)" 
Range("I71").Select 
ActiveCell.FormulaR1C1 = "=SUM(R[-61]C:R[-2]C)" 
Range("172").Select 
Active Window.LargeScroll Down:=-3 
Range("J9"). Select 
ActiveCell.FormulaR1C1 = "Ti*B4" 
Range("J9").Select 
Selection.Font.Bold = True 
With Selection 
-Horizontal Alignment = xlCenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
.Onentation = xlHorizontal 
End With 


ActiveWindow.LargeScroll ToRight:=-1 
Active Window.SmallScroll ToRight:=2 
Range("J10").Select 
ActiveCell.FormulaR1Cl1 = "=RC[-2]*R[-6]C[-8]" 
Range("J10"). Select 
Selection.AutoFill Destination:=Range("J10:J69"), Type:= _ 
xIFillDefault 
Range("J10:J69"). Select 
Selection.Replace What:="“B*", Replacement:="*B4", LookAt:=xlPart, _ 
SearchOrder:=xlByRows, MatchCase:=False 
Range("K9").Select 
Active Window.SmallScroll ToRight:=2 
ActiveCell.FormulaR1Cl1 = "Ti*B6" 
Range("K9").Select 
Selection.Font.Bold = True 
With Selection 
HorizontalAlignment = xJCenter 
.VerticalAlignment = xIBottom 
.WrapTent = False 
.Onentation = xlHorizontal 
End With 
Range("K10").Select 
ActiveCell.FormulaR1C1 = "=RC[-3|4R[~4]C[-9]" 
Range("K10").Select 
Selection. AutoFill Destination:=Range("K10:K69"), Type:= _ 
xlFillDefault 
Range(""K 10:K69"). Select 
Selection.Replace What:="*B*", Replacement:="“B6", LookAt:=xI[Part, _ 
SearchOrder:=xlByRows, MatchCase:=False 
Range("L9"). Select 
ActiveCell.FormulaR1C1 = "Ln TI*Ti*B6" 
Range("L9").Select 
Selection.Font.Bold = True 
With Selection 
-Honzontal Alignment = xlCenterAcrossSelection 
.VerticalAlignment = xlBottom 
.WrapText = True 
Orientation = x]Honzontal 
End With 
Range("L9").Select 
ActiveCell.FormulaR1C1 = "Ln TI * Ti4B6" 
With ActiveCell.Characters(Start:=1, Length:=13).Font 
Name = "Arial" 
FontStyle = "Bold" 
size = 10 
strikethrough = False 
Superscript = False 
Subscript = False 
OutlineFont = False 
Shadow = False 
Underline = xlNone 
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.ColorIndex = x!Automatic 
End With 
Range("L10").Select 
ActiveCell.FormulaR!C1] = "=RC[-3]*RCj-1]" 
Range("L10").Select 
Selection. AutoFill Destination:=Range("L10:L69"), Type:= _ 

xIFillDefault 
Range("L10:L69").Select 
Range("M64").Select 
ActiveWindow.LargeScroll Down:=-3 
ActiveWindow.SmallScroll Down:=12 
Active Window. ScrollRow = 28 
ActiveWindow.ScrollRow = 1 
Active Window. ScrollRow = 61 
Range("K70").Select 
ActiveCell.FormulaR!C1 = "SUM" 
Range("K70").Select 
Selection.Font.Bold = True 
With Selection 

Horizontal Alignment = xlCenter 

.VerticalAlignment = xiBottom 

.WrapText = False 

Orientation = xlHorizontal 
End With 
Selection.Copy 
Range("L70").Select 
ActiveSheet.Paste 
Range("K71").Select 
Application.CutCopyMode = False 
ActiveCell.FormulaR1C1 = "=SUM(R[-61]C:R[-2]C)" 
Range("L71").Select 
ActiveCell.FormulaR1C1 = "=SUM(R[-61]C:R[-2]C)" 
Range("L72").Select 
With Toolbars(8) 

Left = 57 

.Top = 106 
End With 
ActiveWindow.LargeScroll Down:=-4 
Range("M9").Select 
ActiveCell.FormulaR1C1 = "Ti4Q6" 
Range("M9"). Select 
Selection.Font.Bold = True 
With Selection 

Horizontal Alignment = xlCenter 

. VerticalAlignment = x1Bottom 

.WrapText = False 

-Onentation = xlHorizontal 
End With 
Range("M10").Select 
ActiveCell.FormulaR1C1 = "=RC[-5]4R[+]C[4]" 
Range("M10").Select 
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Selection. AutoFill Destination:=Range("M10:M69"), Type:= _ 
xIFillDefault 
Range("M10:M69"). Select 
Selection.Replace What:="“Q*", Replacement:="Q6", LookAt:=xlPart, _ 
SearchOrder:=xlByRows, MatchCase:=False 
ActiveWindow.SmallScroll Down:=8 
Active Window. ScrollRow = 33 
Active Window.ScrollRow = 43 
ActiveWindow.ScrollRow = 48 
ActiveWindow.ScrollRow = 51 
Active Window.ScrollRow = 57 
Range("M/70"). Select 
ActiveCell.FormulaR1C1 = "SUM" 
Range("M70"). Select 
Selection.Font.Bold = True 
With Selection 
HorizontalAlignment = xlCenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
Range("M71").Select 
ActiveCell.FormulaR1C1 = "=SUM(R[-61]C:R[-2]C)" 
Range("M72").Select 
ActiveCell.FormulaR1C1 = "mean" 
Range("M73"). Select 
ActiveCell.FormulaR1C1 = "" 
Range("M72").Select 
Selection.Font.Bold = True 
With Selection 
HorizontalAlignment = xiCenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
Range("M73").Select 
ActiveCell.FormulaR1C1 = "=R[-2}C/60" 
Range("M74").Select 
Active Window. ScrollRow = 1 
Active Window.SmallScroll ToRight:=2 
Range("N9"). Select 
ActiveCell.FormulaR1C1 = "(Ti4Q6-mean)*2" 
Range("N9"). Select 
With Selection 
HorizontalAlignment = xlCenterAcrossSelection 
.VerticalAlignment = xlBottom 
.WrapText = True 
Orientation = xlHorizontal 
End With 
Selection. Font.Bold = True 
ActiveWindow.SmallScroll Down:=59 
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Active Window.LargeScroll Down:=-+ 
Range("N10").Select 
ActiveCell.FormulaR1Cl = "=(RC[-1]-R[63}C[-1])*2" 
Range("N10").Select 
Selection. AutoFill Destination:=Range(""N10:N69"), Type:= _ 
xlFillDefault 
Range("N10:N69").Select 
Selection.Replace What:="-M*)*2", Replacement:="-M73)*2", LookAt _ 
=xlPart, SearchOrder:=xlIByRows, MatchCase:=False 
Range("O9"). Select 
ActiveCell.FormulaR1C1 = "(Ln Ti)*2" 
Range("O9").Select 
Selection.Font.Bold = True 
With Selection 
Horizontal Alignment = xICenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
Range("O10").Select 
ActiveCell. FormulaRIC] = "=RC[-6]*2" 
Range("O10").Select 
Selection. AutoFill Destination:=Range(""O10:069"), Type:= _ 
xlFillDefault 
Range("O10:069").Select 
Range("M70"). Select 
Selection.Copy 
Range("N70").Select 
ActiveSheet.Paste 
Range(""O70").Select 
ActiveSheet.Paste 
Range("N71").Select 
Application.CutCopyMode = False 
ActiveCell.FormulaR1C1 = "=SUM(R[-61jC:R[-2]Q)" 
Range(''O71").Select 
ActiveCell.FormulaR1C1 = "=SUM(R[-61]C:Rj-2]C)" 
Range("0O72").Select 
Active Window. SmallScroll ToRight:=2 
Active Window.LargeScroll Down:=-3 
ActiveWindow. ScrollRow = | 
ActiveWindow.LargeScroll ToRight:=-1 
With Toolbars(8) 
Left = 543 
.Top = 37 
End With 
Range("B4").Select 
ActiveWindow.SmallScroll ToRight:=6 
Active Window. SmallScroll Down:=60 
ActiveWindow.LargeScroll Down:=4 
ActiveWindow.LargeScroll ToRight:=-1! 
Range('"B4").Select 
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ActiveCell.FormulaR1Cl = _ 
"=((6/P1IN*2))*(R[67]C[ 13 ]-(R[67]C[7]*2)/60))*(1/59))*(-0.5)" 
Range("B4"). Select 
Selection.Copy 
ActiveWindow.SmallScroll ToRight:=10 
Range("Q4"). Select 
ActiveSheet.Paste 
Application.CutCopyMode = False 
ActiveCell. FormulaR1Cl = _ 
"=((6/(PIQ*2))*(RI67]C[-2]-(R[67] C[-8 ]42)/60))*(1/59))*(-0.5)" 
Range("Q5").Select 
ActiveWindow.LargeScroll ToRight:=-1 
ActiveWindow.SmallScroll ToRight:=-1 
Range("C6").Select 
ActiveCell.FormulaRICl = _ 
"=(60/RC[-1])+R[65]C[6]-(60*R[65]C[9]/R[65]C[8])" 
Range("C7").Select 
ActiveWindow.LargeScroll ToRight:=-1 
Range("B9").Select 
ActiveCell.FormulaR1 C1 = "=(60/R[62]C[9])4(1/R[-3]©Q)" 
Range("B10").Select 
Active Window.SmallScroll ToRight:=-1 
Range("B11").Select 
ActiveCell.FormulaR1C1 = "=1/R[-2]C" 
Range("B12").Select 
Active Window.SmallScroll ToRight:=1 1 
Range("P6").Select 
ActiveCell. FormulaR1C1 = "=(R[65]C[-3])/((R[65]C[-2])(0.5))" 
Range("Q9"). Select 
Active Window.LargeScroll ToRight:=-2 
ActiveWindow.SmallScroll ToRight:=9 
Range("Q10").Select 
ActiveCell.FormulaR1Cl1 = "theta" 
Range("Q10").Select 
Selection.Font.Bold = True 
With Selection 
HorizontalAlignment = x]Center 
. Vertical Alignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
Range("Q9"). Select 
ActiveCell.FormulaR1C1! = "=1/(((1/60)*R[62]C[-4](1/RI-3]C))" 
Range("Q11").Select 
ActiveCell. FormulaR!I C1 = "=1/R[-2]C" 
Range("Q12").Select 
End With 


End Sub 
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APPENDIX C. MACRO 100 ENGINES 


"MACRO TO 100 ENGINES 
‘ thesis Macro 
‘Macro recorded by FERNANDO GUIMIL 


Sub thesis() 
Sheets. Add 
Sheets("Sheet2").Select 
Range("A1:C2").Select 
ActiveCell.FormulaR1C1 = "MAXIMUM LIKELIHOOD ESTIMATOR" 
Range("A1:C2").Select 
Range("A2"). Activate 
Selection.Font.Bold = True 
With Selection 
-HonizontalAlignment = xlCenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
.Onentation = xlHorizontal 
End With 
With Selection 
-HorizontalAlignment = xlCenterAcrossSelection 
.VerticalAlignment = xIBottom 
.WrapText = True 
.Onentation = xlHorizontal 
End With 
Range("D1").Select 
ActiveWindow.SmallScroll ToRight:=10 
With Toolbars(8) 
Left = 241 
.Top = 228 
End With 
ActiveWindow.LargeScroll ToRight:=-1 
Range("P1:R1").Select 
ActiveCell.FormulaR1C1 = "MODIFIED MOMENT METHOD" 
Range("P1:R1").Select 
Range("Q1"). Activate 
Selection.Font.Bold = True 
With Selection 
-HonzontalAlignment = x]Center 
. VerticalAlignment = xIBottom 
.WrapText = False 
.Onentation = xlHorizontal 
End With 
With Selection 
-HorizontalAhgnment = x!CenterAcrossSelection 
.VerticalAlignment = xlBottom 
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.WrapText = True 
Orientation = x!Honzontal 
End With 
Range("N1").Select 
ActiveWindow.LargeScroll ToRight:=-1 
With Toolbars(8) 
Left = 436 
.Top = 226 
End With 
Range("A2").Select 
ActiveCell.FormulaR1C1 = "shape parameter" 
With ActiveCell.Characters(Start:=1, Length:=15).Font 
.Name = "Arial" 
.FontStyle = "Bold" 
size = 10 
Strikethrough = False 
Superscript = False 
Subscript = False 
-OutlineFont = False 
Shadow = False 
.Underline = xINone 
.ColorIndex = xl Automatic 
End With 
Range("B3").Select 
ActiveCell.FormulaR1C1 = "initial beta" 
Range("B3").Select 
Selection.Font.Bold = True 
With Selection 
-HorizontalAlignment = xICenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
.Onentation = xlHorizontal 
End With 
Range bs: )isclect 
ActiveCell.FormulaR1C1 = "final beta" 
Range("B5").Select 
Selection.Font.Bold = Tre 
With Selection 
.HorizontalAlignment = xICenter 
.VerticalAlignment = xIBottom 
.WrapText = False 
.Onentation = xlHorizontal 
End With 
Range("A7:C7").Select 
ActiveCell.FormulaRIC1 = "scale parameter" 
Range("A7:C7").Select 
Range("B7"). Activate 
Selection. Font.Bold = Tne 
With Selection 
.HorizontalAlignment = xlCenter 
. VerticalAlignment = xlBottom 
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.WrapTent = False 

Orientation = xlHorizontal 
End With 
Range("A7").Select 
ActiveCell.FormulaRIC] = "" 
Range("B7").Select 
ActiveCell.FormulaR!I C1 = “scale parameter" 
Range("B8").Select 
ActiveCell.FormulaRIC1I = "Janda" 
Range("B8").Select 
Selection.Font.Bold = True 
With Selection 

-Horizontal Alignment = xICenter 

.VerticalAlignment = xlBottom 

.WrapText = False 

Orientation = xlHorizontal 
End With 
ActiveWindow.SmallScroll ToRight:=0 
Range("A2:C2").Select 
Selection.Copy 
ActiveWindow.SmallScroll ToRight:=1 1 
Range("Q2").Select 
ActiveSheet.Paste 
Application.CutCopyMode = False 
Selection. Clear 
ActiveWindow.LargeScroll ToRight:=-2 
Range("A2:C2").Select 
Selection. Copy 
ActiveWindow.SmallScroll ToRight:=9 
Range("P2").Select 
ActiveSheet.Paste 
Range("O4").Select 
ActiveWindow.LargeScroll ToRight:=-] 
Range("B3").Select 
Application.CutCopyMode = False 
Selection.Copy 
ActiveWindow.SmallScroll ToRight:=11 
Range("Q3").Select 
ActiveSheet.Paste 
ActiveWindow.LargeScroll ToRight:=-2 
Range("BS5S").Select 
Application.CutCopyMode = False 
Selection. Copy 
ActiveWindow.SmaliScroll ToRight:=12 
Range("Q5").Select 
Application.CutCopyMode = False 
Selection. Copy 
ActiveSheet.Paste 
Range("R7").Select 
ActiveWindow.LargeScroll ToRight:=-2 
Range("B5").Select 
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Application.CutCopyMode = False 
Selection. Copy 
Active Window.SmallScroll ToRight:=12 
Range("Q5").Select 
ActiveSheet.Paste 
Active Window.LargeScroll ToRight:=-2 
Range("A7:C8"). Select 
Application. CutCopyMode = False 
Selection.Copy 
ActiveWindow.SmallScroll ToRight:=1 1 
Range("P7").Select 
ActiveSheet.Paste 
Range("N7").Select 
ActiveWindow.LargeScroll ToRight:=-2 
Range("D11").Select 
Application.CutCopyMode = False 
Selection. Copy 
Range("D7").Select 
ActiveSheet.Paste 
Range("B8").Select 
Active Window.SmallScroll ToRight:=13 
ActiveWindow.LargeScroll ToRight:=-2 
ActiveWindow.SmallScroll ToRight:=14 
ActiveWindow.LargeScroll ToRight:=-2 
Range("E1").Select 
Application.CutCopyMode = False 
ActiveCell.FormulaR1C] = "Number of observations" 
Range("E1").Select 
Selection. Font.Bold = True 
With Selection 
HorizontalAlignment = x!Center 
.VerticalAlignment = xIBottom 
.WrapText = False 
Orientation = xl!Honizontal 
End With 
With Selection 
-HorizontalAlignment = xlCenterAcrossSelection 
. Vertical Alignment = xIBottom 
.WrapText = True 
Orientation = xlHorizontal 
End With 
Range("E2"). Select 
ActiveCell.FormulaR1Cl1 = "n" 
Range("E2").Select 
Selection. Font. Bold = True 
With Selection 
HorizontalAlignment = xiCenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHonzontal 
End With 
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With Toolbars(8) 

Left = 539 

.Top = 30 
End With 
Sheets("Sheet1").Select 
Range("B2:B101").Select 
Selection.Copy 
Sheets("Sheet2").Select 
Active Window.SmallScroll Down:=1 
Range("D 10"). Select 
ActiveSheet.Paste 
Sheets("Sheet1"').Select 
Range("C2:C101").Select 
Application.CutCopyMode = False 
Selection.Copy 
Sheets("Sheet2").Select 
Range("E10").Select 
ActiveSheet.Paste 
Sheets("Sheet 1"). Select 
Range("D2:D101").Select 
Application. CutCopyMode = False 
Selection.Copy 
Sheets("Sheet2").Select 
Range("F10").Select 
ActiveSheet.Paste 
Active Window.LargeScroll ToRight:=0 
Range("D10:D 109"). Select 
Application.CutCopyMode = False 
Selection.Copy 
Range("G10").Select 
ActiveSheet.Paste 
Range("E10:E109").Select 
Application.CutCopyMode = False 
Selection.Copy 
Active Window.SmallScroll Down:=100 
Range("G110").Select 
ActiveSheet.Paste 
Active Window.LargeScroll Down:=-6 
Range("F10:F109"). Select 
Application.CutCopyMode = False 
Selection.Copy 
Active Window.SmallScroll Down:=200 
Range("G210"). Select 
ActiveSheet.Paste 
Active Window.LargeScroll Down:=-11 
Range("G10:G309") Select 
Application.CutCopyMode = False 
Selection. Copy 
Range("H10").Select 
ActiveSheet.Paste 
Application.CutCopyMode = False 
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Selection.Sort Keyl :=Range("H10"), Orderl:=xlAscending, Header:= _ 
xlGuess, OrderCustom:=1, MatchCase:=False, Orientation:= _ 
xlTopToBottom 

ActiveWindow.SmallScroll Down:=-1 

Range("H2:H6"). Select 

ActiveCell.FormulaR1C1 = "failures times Ti" 

Range("H2:H6").Select 

Range("H3"). Activate 

Selection.Font.Bold = True 

With Selection 
.HorizontalAlignment = xlCenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
-Onentation = xlHorizontal 

End With 

With Selection 
HorizontalAlignment = xlCenterAcrossSelection 
.VerticalAlignment = xlBottom 
.WrapText = True 
Orientation = x!Horizontal 

End With 

Range("H9"). Select 

ActiveCell.FormulaR1C1 = "Ti" 

Range("H9").Select 

Selection.Font.Bold = True 

With Selection 
Horizontal Alignment = xICenter 
.VerticalAlignment = xIBottom 
.WrapText = False 
.Onentation = xlHorizontal 

End With 

Range("I9").Select 

ActiveCell.FormulaR1C1 = "Ln Ti" 

Range("I9").Select 

Selection.Font.Bold = True 

With Selection 
-HorizontalAlignment = xlCenter 
.VerticalAlignment = x1Bottom 
.WrapText = False 
-Onentation = xlHorizontal 

End With 

Range("I10"). Select 

ActiveCell. FormulaR1C1 = "=LN(RC[-1])" 

Range("110").Select 

Selection.AutoFill Destination:=Range("110:1309"), Type:= _ 
xIFillDefault 

Range("T10:1309").Select 

With Toolbars(8) 

.Lefi = 79 
.Top = 270 
End With 
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Range("G315").Select 
Active Window.SmallScroll ToRight:=2 
Range("I305").Select 
ActiveWindow.LargeScroll Down:=-17 
Range("J9"). Select 
ActiveWindow.LargeScroll ToRight:=-1 
ActiveWindow.SmallScroll ToRight:=2 
Range("J9").Select 
ActiveCell.FormulaR1C1 = "Ti‘B4" 
Range("J9").Select 
Selection.Font.Bold = True 
With Selection 
-HorizontalAlignment = xlCenter 
.VerticalAlignment = xIBottom 
.WrapTent = False 
Orientation = xlHorizontal 
End With 
Selection.Copy 
Range("K9").Select 
ActiveSheet.Paste 
Active Window.LargeScroll ToRight:=-1 
ActiveWindow.SmallScroll ToRight:=4 
Application.CutCopyMode = False 
ActiveCell.FormulaRIC1 = "Ti*B6" 
Range("J10").Select 
ActiveCell.FormulaR1IC1 = "=(RC[-2]/R{-6]C[-8})" 
Range("J10").Select 
Selection. AutoFill Destination:=Range("J10:J309"), Type:= _ 
xiFillDefault 
Range("J10:J309"). Select 
Range("J300"). Select 
ActiveWindow.LargeScroll Down:=-16 
Active Window.SmallScroll Down:=-9 
Range("J10:J309"). Select 
Selection.Replace What:="“B*)", Replacement:="“B4)", LookAt:= _ 
xlPart, SearchOrder:=x!ByRows, MatchCase:=False 
Range("K 10"). Select 
ActiveCell.FormulaR1Cl1 = "=(RC[-3}4R[~4]C[-9]})" 
Range("K10").Select 
Selection. AutoFill Destination:=Range("K 10:K309"), Type:= _ 
x1FillDefault 
Range("K10:K309"). Select 
Selection.Replace What:="“B*)", Replacement:="*B6)", LookAt:= _ 
xlPart, SearchOrder:=xIByRows, MatchCase:=False 
Range("L8:L9"). Select 
ActiveCell.FormulaR1C1 = "Ln Ti * Ti*B6" 
Range("L8:L9"). Select 
Range("L9"). Activate 
With Selection 
-HorizontalAlignment = xlCenterAcrossSelection 
. VerticalAlignment = xlBottom 


se) 


.WrapText = True 
Orientation = xlHorizontal 
End With 
Range("L9").Select 
Selection.Delete Shift:=xlUp 
Range("L8").Select 
Selection. Copy 
Range("L9").Select 
ActiveSheet. Paste 
Range("L8").Select 
Application. CutCopyMode = False 
Selection.Clear 
Range("L9").Select 
Selection.Font.Bold = True 
With Selection 
HorizontalAlignment = xlCenter 
.VerticalAlignment = xIBottom 
.WrapText = True 
Orientation = xlHorizontal 
End With 
Range("L10").Select 
ActiveCell.FormulaR1C1 = "=RC[-3]*RC[-1]" 
Range("L10").Select 
Selection. AutoFill Destination:=Range("L10:L309"), Type:= _ 
xIFillDefault 
Range("L10:L309"),. Select 
Active Window. LargeScroll Down:=-1 
Range("M285").Select 
Active Window. LargeScroll Down:=-16 
ActiveWindow.SmallScroll ToRight:=4 
ActiveWindow.LargeScroll ToRight:=-1 
Active Window.SmallScroll ToRight:=8 
Active Window.SmallScroll Down:=33 
Active Window.LargeScroll Down:=-2 
Active Window.ScrollRow = 174 
ActiveWindow.ScrollRow = 218 
Active Window.ScrollRow = 252 
Active Window. ScrollRow = 309 
Range("L310").Select 
ActtveCell.FormulaR1C1 = "SUM" 
Range("L3 10").Select 
Selection.Font.Bold = True 
With Selection 
-HorizontalAlignment = xlCenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
Range("L311").Select 
ActiveCell. FormulaR1C1] = "=SUM(R[-301]C:R[-2]O)" 
Range("L312").Select 
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Range("M9").Select 
ActiveCell.FormulaR1C1 = "Ti*Q6" 
Range("M9"). Select 
Selection.Font.Bold = True 
With Selection 
-HonizontalAlignment = xICenter 
.VerticalAlignment = xiBottom 
.WrapTent = False 
.Onentation = xlHonizontal 
End With 
Range("M10").Select 
ActiveCell.FormulaR1C1 = "=(RC[-5]4R{-4]C[4])" 
Range("M10").Select 
Selection. AutoFill Destination:=Range("M10:M309"), Type:= _ 
xlFillDefault 
Range("M10:M309"). Select 
Selection.Replace What:="“Q*)", Replacement:="“Q6)", LookAt:= _ 
xlPart, SearchOrder:=xlByRows, MatchCase:=False 
Range("N9").Select 
ActiveWindow.SmallScroll Down:=107 
ActiveWindow.ScrollRow = 162 
ActiveWindow.ScrollRow = 214 
ActiveWindow.ScrollRow = 243 
ActiveWindow.ScrollRow = 261 
Active Window. ScrollRow = 282 
ActiveWindow.ScrollRow = 297 
Range("M310").Select 
ActiveCell.FormulaR1C1 = "SUM" 
Range("M310").Select 
Selection.Font.Bold = True 
With Selection 
HonzontalAlignment = xICenter 
.VerticalAlignment = xIBottom 
.WrapText = False 
.Onientation = xlHorizontal 
End With 
Range("M311").Select 
ActiveCell.FormulaR1C1 = "=SUM(R[-301]C:R[-2]©)" 
Range("M3 12").Select 
ActiveCell.FormulaR1C1 = "mean" 
Range("M3 12").Select 
Selection.Font.Bold = True 
With Selection 
-HonzontalAlignment = x!Center 
.VerticalAlignment = xlBottom 
.WrapText = False 
.Onientation = xlHonzontal 
End With 
Range("M313").Select 
ActiveWindow.LargeScroll Down:=-17 
Range("H3").Select 
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Active Window.LargeScroll ToRight:=-|! 
Range("E3").Select 
ActiveCell.FormulaR1C1 = "300" 
Range("E3").Select 
Selection.Font.Bold = True 
With Selection 
Horizontal Alignment = xiCenter 
.VerticalAlignment = xIBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
ActiveWindow.SmallScroll ToRight:=5 
ActiveWindow.ScrollRow = 144 
Active Window. ScrollRow = 178 
ActiveWindow.ScrollRow = 219 
Active Window. ScrollRow = 264 
ActiveWindow.ScrollRow = 291 
ActiveWindow.ScroliRow = 312 
ActiveWindow.ScrollRow = 274 
ActiveWindow.LargeScroll Down:=-15 
Active Window.SmallScroli Down:=37 
ActiveWindow.ScrollRow = 283 
Active Window.ScrollRow = 312 
Active Window. ScrollRow = 303 
Range("M313").Select 
ActiveCeli.FormulaR1C1 = "=R[-2]C/300" 
Range("M314").Select 
ActiveWindow.LargeScroll Down:=-17 
Range("N9"). Select 
ActiveCell.FormulaR1C1 = "(Ti*Q6-mean)*2" 
Range("N9").Select 
Selection.Font.Bold = True 
With Selection 
.HorizontalAlignment = xlCenter 
. VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
With Selection 
Horizontal Alignment = xlCenterAcrossSelection 
.VerticalAlignment = xlBottom 
.WrapText = True 
Orientation = xiHorizontal 
End With 
Range("N10"). Select 
ActiveCell.FormulaRIC1 = "=(RCJ-2]-R[303]C[-1])" 
Range("N10"). Select 
Selection. AutoFill Destination:=Range("N10:N309"), Type:= _ 
xIFillDefault 
Range("N10:N309"). Select 
ActiveWindow.SmallScroll ToRight:=1 
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Range("N10").Select 
ActiveCell.FormulaR1C1 = "=(RC]-1]-R[303]C{[-1])*2" 
Range("N10").Select 
Selection. AutoFill Destination:=Range("N10:N309"), Type:= _ 
xlFillDefault 
Range("N10:N309").Select 
Selection.Replace What:="-M*)*2", Replacement:="-M313)*2", LookAt _ 
=x]Part, SearchOrder:=xIByRows, MatchCase:=False 
Active Window.SmallScroll Down:=-2 
ActiveWindow.SmallScroll ToRight:=3 
ActiveWindow.LargeScroll ToRight:=-1 
Range("B4").Select 
ActiveCell.FormulaR1C1 = "=((6/(PIQ%2))*1)" 
Range("B5").Select 
Active Window. SmallScroll Down:=16 
Active Window. ScrollRow = 1 
Active Window.SmallScroll Down:=20 
ActiveWindow.ScrollRow = 612 
Active Window. ScrollRow = 581 
ActiveWindow.ScrollRow = 539 
ActiveWindow. ScrollRow = 452 
ActiveWindow.LargeScroll Down:=-9 
Active Window. ScrollRow = 339 
ActiveWindow. ScrollRow = 309 
Range("H310").Select 
ActiveCell.FormulaR1C1 = "SUM" 
Range(""H3 10").Select 
Selection.Font.Bold = True 
With Selection 
.HorizontalAlignment = xlCenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
Range("H311").Select 
ActiveCell.FormulaR1C] = "=SUM(R[-301]C:R[-2]C)" 
Range("H312").Select 
ActiveWindow. ScrollRow = 1 
ActiveWindow. ScrollRow = 465 
ActiveWindow. ScrollRow = 287 
ActiveWindow. ScrollRow = 322 
Active Window.LargeScroll ToRight:=1 
Active Window. ScrollRow = 203 
ActiveWindow.ScrollRow = 269 
Active Window.ScrollRow = 3 
Active Window.SmallScroll ToRight:=-4 
ActiveWindow.ScroliRow = 124 
ActiveWindow. ScrollRow = 223 
Active Window. ScrollRow = 317 
Active Window. ScrollRow = 274 
Active Window. ScrollRow = 276 
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ActiveWindow.SmallScroll Down:=23 
Range("H310").Select 
Selection.Copy 
Range("I3 10").Select 
ActiveSheet. Paste 
Range("I311").Select 
Application.CutCopyMode = False 
ActiveCell.FormulaR1C1 = "=SUM(R[-301]C:R[-2]C)" 
Range("I312").Select 
Active Window.ScrollRow = 1 
ActiveWindow.SmallScroll Down:=1 
ActiveWindow.SmallScroll ToRight:=2 
Range("O9").Select 
ActiveCell.FormulaR!C1 = "(Ln Ti)%2" 
Range("O9").Select 
Selection.Font.Bold = True 
With Selection 
-HorizontalAlignment = xICenter 
. VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
Range("O10").Select 
ActiveCell.FormulaR1C1 = "=RC[-5]*2" 
Range("O10").Select 
Selection. AutoFill Destination:=Range("O10:0309"), Type:= _ 
xlFillDefault 
Range("O10:0309"). Select 
Range("O310").Select 
ActiveCell.FormulaR1C1 = "SUM" 
Range("O310").Select 
Selection.Font.Bold = True 
With Selection 
Horizontal Alignment = xICenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
Range("O311").Select 
ActiveCell.FormulaR1C1 = "=SUM(R[-301]C:R[-2]C)" 
Range("O312").Select 
Active Window.SmallScroll ToRight:=-7 
Active Window.LargeScroll Down:=-15 
ActiveWindow.ScrollRow = | 
Active Window.SmallScroll ToRight:=1 
Range("B4"). Select 
ActiveCell FormulaR1C] = _ 
"=((6/(P10%2))*(R[307]C[13]-R[307]C[7]/300)/299)(-0.5)" 
Range("B4").Select 
ActiveWindow.SmallScroll ToRight:=6 
ActiveWindow.SmallScroll Down:=- 1 
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ActiveCell.FormulaR1IC1 = _ 
"=((6/(PIQ%2))*(R[307]C| 13 ]-(R[307]C[7]*2)/300V/299)(-0.5)" 

Range("B+4").Select 

Active Window. ScrollRow = 419 

ActiveWindow.ScrollRow = 300 

ActiveWindow.SmallScroll ToRight:=7 

ActiveWindow.ScrollRow = 1 

Active Window.LargeScroll ToRight:=- 1 

ActiveCell.FormulaR1Cl = _ 
"=((6/(PIQ%2))*(R[307]C[13]-(R[307]C[7]*2)/300)/299)(-0.5)" 

Range("B4").Select 

ActiveCell.FormulaR1Cl] = _ 
"=((6/(PIQ%2))*(R[307] C[13 ]-C(R[307] C7 }42)/300)/299)(-0.5)" 

Range("I10").Select 

ActiveWindow.SmallScroll ToRight:=2 

Range("O10").Select 

Active Window.LargeScroll ToRight:=-1 

Range("B4").Select 

ActiveCell.FormulaR1Cl = _ 
"=((6/(PIQ%2))*(R[307]C[13]-CR[307]C[7]*2)/300)/299)(-0.5)" 

Range("B+").Select 

ActiveWindow.SmallScroll ToRight:=2 

Range("J10").Select 

Active Window.LargeScroll ToRight:=-1 

ActiveWindow.SmallScroll ToRight:=2 

Active Window.LargeScroll ToRight:=-1 

Range("B4").Select 

ActiveWindow.SmallScroll ToRight:=2 

Range("J11").Select 

ActiveWindow.LargeScroll ToRight:=-1 

Range("B4"). Select 

ActiveCell.FormulaR1Cl = _ 
"=((6/(PIQ%2))*(R[307] C[13]-CR[307]C[7}42)/300)/299)(-0.5)" 

Range("O10").Select 

ActiveWindow.LargeScroll ToRight:=-1 

ActiveWindow.SmallScroll ToRight:=7 

ActiveCell.FormulaR1C1 = "=RC[-6]42" 

Range("O10"). Select 

Selection. AutoFill Destination:=Range("O10:0309"), Type:= _ 
xIFillDefault 

Range("O10:0309").Select 

ActiveWindow.LargeScroll Down:=-16 

ActiveWindow.SmallScroll Down:=-12 

ActiveWindow.LargeScroll ToRight:=-] 

Range("B4").Select 

ActiveWindow.SmallScroll ToRight:=3 

Active Window.LargeScroll ToRight:=-1 

Active Window.SmallScroll ToRight:=13 

Range("Q4"). Select 

ActiveCell.FormulaR1C1 = "“=RC[-15]" 

Range("Q5").Select 


6] 


ActiveWindow.LargeScroll ToRight:=-2 
With Toolbars(8) 
Lent = 39% 
.Top = 139 
End With 
ActiveWindow.LargeScroll ToRight:=1 
ActiveWindow.SmallScroll ToRight:=-9 
ActiveWindow.SmallScroll Down:=3 
Range("C6").Select 
ActiveCell.FormulaR1Cl1 = _ 
"=(300/RC[-1])+R[305]C[6]-(300*R[305]C[9/R[305]C[8])" 
Range("K3 10").Select 
ActiveCell.FormulaR1C1 = "SUM" 
Range("K310").Select 
Selection.Font.Bold = True 
With Selection 
HorizontalAlignment = xICenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
Onientation = xlHorizontal 
End With 
Range("K311").Select 
ActiveCell.FormulaR1 C1] = "=SUM(RJ-301)C:R[-2]C)" 
Range("K3 12").Select 
ActiveWindow.LargeScroll Down:=-17 
ActiveWindow.LargeScroll ToRight:=-1 
Range("C6").Select 
ActiveCell.FormulaR1C1 = _ 
"=(300/RC[-1])+R[305]C[6]-(300*R[305]C]9J/RI305]C[8])" 
Range("C6").Select 
ActiveWindow.SmallScroll ToRight:=2 
Active Window. ScrollRow = 432 
ActiveWindow.ScrollRow = 311 
Active Window. ScrollRow = | 
ActiveWindow.SmallScroll ToRight:=-2 
Range("B10").Select 
ActiveCell.FormulaR1C1 = "theta" 
Range("B10").Select 
Selection. Font.Bold = True 
With Selection 
.HorizontalAlignment = x!Center 
.VerticalAlignment = xIBottom 
.WrapText = False 
.Onentation = xlHorizontal 
End With 
Range("B9").Select 
ActiveCell.FormulaR1ICl = "=(300/R[302|C[9])*C1/RI-3]C)" 
Range("B9").Select 
ActiveCell.FormulaR1C1 ="" 
Range("C9").Select 
ActiveCell.FormulaR1Cl = "" 


Range("B9"). Select 
ActiveCell.FormulaR1C1 = "=(300/R[302}C[9])*U/RI-3]C)" 
Range("B11").Select 
ActiveCell.FormulaR1C1] = "=1/R[-2]C" 
Range("B11").Select 
ActiveWindow.SmallScroll ToRight:=9 
Range("Q6").Select 
With Toolbars(8) 
Left = 58 
Top =25 
End With 
Range("P6"). Select 
ActiveCell.FormulaR1C1 = "((1/299)" 
Range("S319").Select 
Active Window. ScrollRow = 302 
Range("N310").Select 
ActiveCell.FormulaR!Cl = "SUM" 
Range("N3 10"). Select 
Selection.Font.Bold = True 
With Selection 
Horizontal Alignment = xlCenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
Range("N311").Select 
ActiveCell.FormulaR1C1 = "=SUM(R[-301]C:R[-2]C)" 
Range("N312").Select 
ActiveWindow.LargeScroll Down:=-17 
Range("P6").Select 
ActiveCell.FormulaR1C1 = "(M311)/((N311)4(0.5))" 
Range("P6").Select 
ActiveWindow.LargeScroll ToRight:=0 
ActiveCell.FormulaR1C1 = _ 
"=R[305]C[-3 })/(R[305]C[-2])*(0.5)" 
Range("Q9").Select 
ActiveCell. FormulaR1 C1 = "=1/(((1/300)*R[302]C[-4])C1/R[-3]C))" 
Range("Q10").Select 
ActiveWindow.LargeScroll ToRight:=-1 
ActiveWindow.SmallScroll ToRight:=13 
ActiveCell.FormulaR!C1 = "theta" 
Range("Q10").Select 
Selection.Font.Bold = True 
With Selection 
-HorizontalAlignment = xICenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End With 
Range("Q11").Select 
ActiveCell.FormulaR1C1 = "=1/R[-2]C" 
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Range("Q12").Select 
ActiveWindow.LargeScroll ToRight:=-2 
End Sub 
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APPENDIX D. PARAMETER CALCULATIONS FOR THE 20 AND 100 ENGINES 
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APPENDIX E. FUNCTION “tocompare” 


function(v, betamle, thetamle, betammm, thetammm) 
{ # This function calculates the kolmogorov-smimov goodness of the fit, 
# and plots of the empical distribution and theoretical distribution 
# fora set of data and two sets of parameters of weibull distribution. 
win. graphQ) 
a <- ks.gof(v, dist = "weibull", shape = betamle, scale = thetamle) 
b <- ks.gof(v, dist = "weibull", shape = betammm, scale = thetammm) 
par(mfrow = c(1, 2)) 
cdf.compare(v, dist = "weibull", shape = betamle, scale = thetamle) 
cdf.compare(v, dist = "weibull", shape = betammm, scale = thetammm) 
list(a, b) 
} 


12 


tocompare(a20. 1,0.973,251.71,0.9857,253. 14) 
[U1]: 


One-sample Kolmogorov-Smimov Test; hypothesized distribution = weibull 


data: v 
ks = 0.0546, p-value = 0.994 


alternative hypothesis: True cdf does not equal the weibull Distn. for at least one sample point. 


([2]]: 


One-sample Kolmogorov-Smimov Test; hypothesized distribution = weibull 
data: v 
ks = 0.0509, p-value = 0.9977 


Alternative hypothesis: True cdf does not equal the weibull Distn. for at least one sample point. 
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Figure 5. Kolmogorov-Smirnov Test; Hypothesized Distribution = Weibull. (a) ks = 
0.0546, p-value = 0.994, (b) ks = 0.0509, p-value = 0.9977 
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tocompare(al00.1,1.261,292.56,1.2882,294 66) 
[C1]: 


One-sample Kolmogorov-Smirmov Test; hypothesized distribution = weibull 
data: v 


ks = 0.0319, p-value = 0.9207 


Alternative hypothesis: True cdf does not equal the weibull Distn. for at least one sample point. 


[2]: 


One-sample Kolmogorov-Smimov Test; hypothesized distribution = weibull 


data: v 


ks = 0.0382, p-value = 0.774] 


Alternative hypothesis: True cdf does not equal the weibull Distn. for at least one sample point. 
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Figure 6. Kolmogorov-Smirnov Test; Hypothesized Distribution = Weibull. (a) ks = 
0.0319, p-value = 0.9207, (b) ks = 0.0382, p-value = 0.7741 
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